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ABSTRACT 

We investigate the effect of the hot CNO cycle breakout reaction ls O(a,7) 19 Ne on the occurrence of type I 
X-ray bursts on accreting neutron stars. For < 0. 1 , w here is a dimensionless factor by which we multiply 
the ls O(a,7) 19 Ne reaction rate of Caugh lan & Fowlerl fl988). our model predicts that bursts should occur only 
for accretion rates M below a critical value w 0.3MEdd- This agrees with observations. For larger values of 
including the standard choice f w = 1, the model switches to a new regime in which bursts occur all the way 
up to M ss MEdd- Since the latter regime disagrees with observations, we suggest that the true 15 0(a,7) 19 Ne 
reaction rate is lower than usually assumed. 

Subject headings: dense matter — nuclear reactions, nucleosynthesis, abundances — stars: neutron — X-rays: 
binaries — X-rays: bursts 



1. INTRODUCTION 

Type I X-ray bursts are thermonuclear explosions that 
occur on accreting neutron stars in low-mass X-ray bina- 
ries. They are triggered by thermally unstable H or He 
burning near the stellar surface (for reviews, see ICum ming 
120041 Strohmav er & Bildstenll2006h . For systems with ac- 
cretion rates M < O.lMEdd, where MEdd denotes the mass 
accretion rate at which the accretion luminosity is equal to 
the Eddington limit, the basic physics of the burst onset is 
well understoo d to be that of the thin shell thermal in stabil- 
ity ( Schwarzschild & Harmll965UHansen & van Hornll975l) . 
and theoretical models have been rather successful at re- 
producing the g ross characteristics of burst observations in 
this regime (e.g.. |Fuiimoto et alJl98lllFushi ki & Lambll987t 
ICumming & Bildstenl2000tlNaravan & Hevll2003l) . 

This Letter addresses a longstanding problem that afflicts 
nearly all burst models when the accreti on rate M > O.IMfm- 
Both simp le one-zone burst models ( Fuiimoto et al J ll 98 It 
lPaczvnskaU983l lBildstenlH998l: iHeger et al J 120051) and so- 
phisti cated time- depend ent multi-zone mo dels (|Ay asli & Jossl 
ll982HTaam et al]ll996tlFisker et alJl200ltlHeger et al.l2005l) 
predict that bursts should occur for all M up to « MEdd- 
Observations, however, in dicate that bursts do not occur for 
M > 3Mbhh (Ivan Paradiis et al.ll979l ll988: Cornelisse et al 
120031 iRemi llard et alJ l2006l) . Furthermore. Ivan Parad iis et al 
(119881) found that, for 0.1 <M /M Edd < 0.3, a significant frac- 
tion of the accreted plasma burns stably between consecutive 
bursts, leading to large values > 1000 of the parameter a, the 
accretion energy released between successive bursts divided 
by the nuclear energy released during a burst. Most theoret- 
ical models, on the other hand, predict that nearly all of the 
accreted matter burns unstably during bursts, giving a nearly 
c onstanta < 100 a t all M 

iNaravan & Hevll J2003L hereafter NH03) developed a 
global linear stability analysis of the accreted plasma on the 
surface of a neutron star. They discovered a new regime of 
unstable nuclear burning for 0.1 < At/ M^dA < . 3 that they re- 
ferred to as "delayed mixed bursts." Their model reproduced 
both the occurrence of considerable stable burning preced- 
ing a burst (leading to values of a > 1000) and the absence 
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of bursts for At > 0.3MEdd, in agreement with observations. 
However, the complexity of their model made it difficult to 
understand the basic physics behind delayed mixed bursts and 
to identify the reasons why their results differed from those of 
othe r theoretical models at hig h accretion rates. To remedy 
this, ICooper & Naravanl (120061 hereafter CN06) constructed 
a simple two-zone model that helped elucidate the physics of 
delayed mixed bursts. They showed that the competition be- 
tween nuclear heating due to triple-a reactions and hot CNO 
cycle H burning on the one hand, and radiative cooling via 
photon diffusion and emission on the other hand, drives an 
overstability that eventually triggers a thin-shell thermal in- 
stability and hence a delayed mixed burst. They asserted that 
H burning via the temperature-independent hot CNO cycle, 
augmented by the extra seed nuclei produced from stable He 
burning, significantly lowers the temperature sensitivity of the 
total nuclear energy generation rate to such an extent as to 
suppress the thin-shell thermal instability for M > 0.3M£dd- 
Therefore, H burning via the hot CNO cycle is ultimately re- 
sponsible for the lower critical At above which bursts do not 
occur in nature. 

The above argument perhaps explains why one-zone mod- 
els fail, since the models generally focus only on He burn- 
ing and make large approximations with respect to H burn- 
ing. But why do detailed time-dependent multi-zone mod- 
els with large reaction networks also perform poorly in re- 
lation to observations? These models are much more so- 
phisticated than the models of NH03 and CN06 and con- 
sequently ought to perform better, whereas in fact the lat- 
ter models agree much better with observations. CN06 hy- 
pothesized that the time-dependent multi-zone burst models 
may have used too large a rate for the experimentally poorly- 
constrained JiotCNO_£y_cJe_bre^out reaction ls O(a,7) 19 Ne 
(Wagoner 1969; Wallace & Woosl evl 1198 It | Langanke et alJ 
119861 IWiescher et al.lll999t iFisker et alJl2006l) . We test this 
hypothesis. We begin in <j2]with a description of the model, 
and we present the results of the model in ^3] We discuss the 
results in |4] and we conclude in ^5] 

2. THE MODEL 

We u se the gen eral-relativistic global linear stability anal- 
ysis of Cooper & Naravan (2005), which is an expanded and 
improved version of the model of NH03, to determine the sta- 



2 



bility of nuclear burning on accreting neutron stars. We as- 
sume that matter accretes spherically onto a neutron star of 
gravitational mass M = 1.4M© and areal radius R = 10.4 km 
at a rate M, where M is the rest mass accreted per unit time 
as measured by an observer at infinity. We set the compo- 
sition of the accreted matter to be that of the Sun, such that 
at the neutron star surface the H mass fraction X oul = 0.7, He 
mass fraction y out = 0.28, CNO mass fraction ZcNO.out =0.016, 
and heavy element fraction Z out = 0.004, where Z refers to all 
metals other than CNO. In this section, we describe the mod- 
ifications to the theoretical model. 

For temperatures T > 8 x 10 7 K, H 
burns predominantly via the hot CNO cycle 
12 C(p,7) f3 N(p,7) 14 0(^ + z/) 14 N(/?,7) 15 0(/3 + z/) 15 N(p,a) 12 C, 
the rate of which is determined by the sl ow /3-decays of 
14 0(fi /9 = 70.6 s) and I5 0(fi /2 = 122 s) llHovle & Fowled 
1965). During the hot CNO cycle, essentially all of the 
CNO ions are converted to I4 and ls O. For the range of 
accretion rates we consider in this investigation, the primary 
br eakout reaction from the ho t CNO cycle into th e rp-process 
oflWallace & Wooslevl (119811) is 15 0(a.^ 1 9 Ne (jGorres et al. 
11 9951 lHahn et al JI1 996t ISchatz et al] 1999). If this breakout 
reaction occurs, the l9 Ne(t\/2 = 17.2 s) produced in this 
reaction can /3-decay and return to the hot CNO cycle 
via I9 Ne(/3 + z/) 19 F(/?, a) l6 0(p, ~f) 17 F(p, 7) 18 Ne(/? + i/) 18 F(p, 
a) 15 0, the effect of which is only to expedite the hot CNO 
cycle by a factor < 1.6. We ignore this effect in our model. 
If I9 Ne captures a proton, however, the resulti ng 20 Na ion 
can n ever return to the hot CNO cycle ( Wallace & Wooslev 
1981). The fate of a Ne ion thus depends on whether or 
not it captures a proton before it can /3-decay. Therefore, the 
probability that a 19 Ne ion is removed from the hot CNO 
cycle is 
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where p is the density, X is the H mass fraction, A P7 ( 19 Ne) 
is the 19 Ne(p,7) 20 Na reaction rate from ICaughlan & Fowler! 
(I1988L hereafter CF88) with electron screening from 
iDewitt et all fl973l and A/3 + ( 19 Ne) = ln2/f 1/2 is the 19 Ne (3- 
decay rate. If 19 Ne captures a proton, further proton captures 
will ensue until the nuclear flow reaches the first waitin g point 
of th e rp-process, which we assume to be 24 Si dWiescher et al.l 
1998). The ls O(a,7) 19 Ne reaction is the slowest reaction of 
this sequence, and so its rate g overns the total reactio n rate of 
this flow. Therefore, we follow Wiesch eir et alJ(H9 99) and ap- 
proximate the total nuclear energy generation rate of the hot 
CNO cycle breakout reaction sequence by 
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where E* = grp/24m p = 1.2 x 10 18 ergsg Y is the He mass 
fraction, Zcno is the mass fraction of CNO elements, and 
A Q7 ( ls O) is the 15 0(o:,7) 19 Ne re action rate from CF88 which 
is based on the rate derived by Lanaanke et al. ( 1986). For 
the temperatures T < 6 x 10 8 K we consider in this work, the 
resonant contribution of the 4.03 Me V state in 19 Ne domi - 
nates the 15 0(a,7) I9 Ne reaction rate (Langankeet al. 1986). 
The reaction rate of this contribution is directly proportional 
to the a width T a . However, r Q is very difficult to mea- 
sure experimentally, and so the true A Q7 ( 15 0) is highly un- 
certa in (for recent progress on experimenta l constraints, see, 
e.g., iTan et al 1 l200l iKanungo et aD l2006). To account for 



this uncertainty, we multiply the CF88 rate by a dimension- 
less free parameter /ip such that / lp = 1 corresponds the the 
usual rate of CF88. We omit the electron screening contribu- 
tion to the reaction rate since the enhancement factor itself 
depen ds on the relatively unconstrained resonance strength 
(Salpeter & van Horn 1969; Miller 119771) . Instead, we sim- 
ply absorb this contribution into fm. Equations (6), (8), (9), 
and (10) of Cooper & Naravan (2005) thus become 
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See Cooper & Naravan (2005) for the definitions of the var- 
ious symbols. Note that the factor of 5 in equation © ac- 
counts for the five proton captures needed to reach 24 Si. For 
a substantial portion of the accretion rates we consider in this 
work, the nuclear flow due to stable burning extends well be- 
yond the relatively small num ber of isotopes included in ou r 
limited network (e.g., Schatz et al. 1999; Fiskeret al. 2005). 
However, the contributions from H and He burning dominate 
the total nuclear energy generation rate and set the thermal 
profile of the accreted layer, and we calculate these contri- 
butions to high accuracy. Consequently, the thermal profiles 
resulting from our network compare very well with those of 
ISchatz et alJ (U99 9) for the range of accretion rates in which 
our calculations overlap. Furthermore, we are interested only 
in the physics of the burst onset, for which only hot CNO cy- 
cle, triple-a, and breakout reactions are important. Thus, our 
limited network should be adequate for our purposes. 

3. RESULTS 

Figure n shows the values of the dimensionless quantity a 
as a function of the Eddington-scaled accretion rate Z acc = 
M/Mndd for five different values of f lv , where a is defined 
as the accretion energy released between successive bursts di- 
vided by the nuclear energy released during a burst. We fol- 
low NH03 and call a type I X-ray burst "prompt" if a < 100 
and "delayed" if a ^ 100. The case /rp = omits the break- 
out reaction sequence entirely, and so this calculation is rep- 
resentative of the results of NH03 and CN06. The models 
gives prompt bursts for Z acc < 0.15 and delayed bursts over 
the range Z acc « 0.15-0.3. For low values of /ip < 0.1, the 
effects of breakout reactions from the hot CNO cycle on the 
occurrence of type I X-ray bursts are minor, and the models 
of NH03 and CN06 continue to provide an adequate descrip- 
tion of the nuclear physics that precedes bursts. However, the 
situation changes for ^0.1. The critical / acc above which 
delayed mixed bursts cease now significantly decreases, and 
the range of accretion rates in which delayed mixed bursts 
occur becomes notably truncated as well. More importantly, 
a new burst regime appears at accretion rates close to Mndd, 
where the bursts have a < 100 and are hence prompt. This 
regime of bursts near the Eddington limit was hypothesized 
by CN06, but it was not present in either the NH03 or CN06 
models due to the exclusion of hot CNO cycle breakout re- 
actions. The usual delayed mixed burst regime and this new 
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FIG. 1. — Plot shows a-values of type I X-ray bursts as a function of the 
Eddington-scaled accretion rate 4cc = M/M£dd for five different values of 
/rp, the dimensionless factor by which we multiply the 15 0(a,7)"Ne reac- 
tion rate from CF88. 



prompt mixed burst regime are separated by a short range of 
accretion rates in which bursts do not occur. These results ob- 
tained with larger v alues of f m agree rath er well wi th previous 
theoretical models (Fuiimoto et al. 1981; Avasli & Joss 1982; 
Paczvriskill9l^lTaamll985UTaam et alJl996llBildstenll998l 
Fisker et al. 2003; Heger et al. 2005), most of which predicted 
bursts to occur up to roughly the Eddington limit. Figure 
^thus explains the origin of the differences between those 
models and the models of NH03 and CN06. If the break- 
out reaction rate is large, say /p > 0.1, the results agree with 
most published models (which include breakout reactions in 
full strength), and if the rate is small, / lp < 0.1, the results 
are similar to those obtained by NH03 and CN06 (who ef- 
fectively set /rp = 0). We reiterate that the latter results agree 
much better with observations. 

4. DISCUSSION 

For / acc < 0.1, the temperatures reached during steady-state 
nuclear burning in the accreted layer are too low for signifi- 
cant leakage out of the hot CNO cycle via breakout reactions, 
and so these reactions have a negligible effect on the onset of 
type I X-ray bursts, regardless of the true ls O(a,7) 19 Ne cross 
section. However, breakout reactions should affect the nuclear 
flow during the burst itself, and hence they could influence 
the bur st lightcurve and possibly also the onset of subsequent 
bursts Jl-isker el alj l2006). CN06 showed that the delayed 
mixed bursts of NH03 occur when Sh, the column depth at 
which H is depleted via stable nuclear burning as measured 
from the stellar surface, is less than but close to Shc the col- 
umn depth at which He is depleted via stable nuclear burning. 
The hot CNO cycle breakout reaction sequence (i) eliminates 
seed nuclei from the hot CNO cycle which slows H burning, 
thereby increasing Eh, and (ii) provides an additional path- 
way by which He may burn, thereby decreasing Sne- Con- 
sequently, one expects the regime of delayed mixed bursts to 
occur at lower l acc if breakout reactions are included. Fur- 
thermore, CN06 showed that the interplay between H burning 
via the hot CNO cycle and He burning via triple-a reactions 
is integral to generating the oscillations that precede delayed 
mixed bursts. Breakout reactions diminish this interplay by 
eliminating the hot CNO cycle seed nuclei and should there- 



fore reduce the range of accretion rates over which delayed 
mixed bursts occur. Figure ^ illustrates both of these effects 
(compare = 1 with /p = 0). 

For accretion rates above those at which delayed mixed 
bursts occur, NH03 found that nuclear burning is always sta- 
ble an d therefore bursts do not occur, in accord with observa- 
tions (Ivan Paradiis et ail 1 19791 fT988t ICornelisse et alJl200l 
Remillard et al. 2006). According to their model, for 0.3 < 
4cc S 1' steady-state H burning via the hot CNO cycle in- 
creases the effective radiative cooling rate and thereby sup- 
presses a He-triggered thin-shell thermal instability (CN06). 
However, their model did not include breakout reactions. If 
the 15 0(a,7) 19 Ne breakout reaction rate is significant in this 
regime, this reaction will suppress hot CNO cycle H burning, 
and thus He burning via triple-a reactions will govern the total 
reaction rate, since the triple-a reaction rate will be the slow- 
est rate in the nuclear flow at these high Z acc . The notion that 
the triple-a reaction rate is the slowest rate in the nuclear flow 
is precisely what one-zone type I X-ray burst ignition mod- 
els assume JFuiim oto et alJ|198U [Paczvrisk illl983t iBildsten 
U998tlCumming & Bildstenl2000UHeger et alJ2005l) . and we 
suggest that this is the reason why the results of one-zone and 
time-dependent multi-zone models agree so well. All of these 
models predict that nuclear burning is thermally stable for 
temperatures T > 5 x 10 8 K, which are reac hed only at ac- 
cretion rates Z acc > 1 (e.g., Schatz et al. 1999). FigM-e^iHus- 
trates that, for relatively large values of there is a regime of 
prompt bursts at accretion rates greater than the accretion rates 
at which delayed mixed bursts occur, and the critical accretion 
rate above which these prompt mixed bursts cease is roughly 
MEdd, in very good agreement with these other models. Fig- 
ure^iUustrates also that nuclear burning is stable for the small 
range of Z acc between these two regimes. These accretion rates 
are high enough to suppress delayed mixed bursts, but they are 
not high enough to cause a sufficient leakage out of hot CNO 
cycle H burning and trigger a thermal instability. 

The model we present in this work suggests that, if the true 
15 0(a,7) 19 Ne cross section is greater than approximately 0.1 
of the CF88 rate, type I X-ray bursts should occur in systems 
with accretion rates near the Eddington limit. Observations 
indicate that this is not the case, since low-mass X-ray bi- 
naries with l acc > 0.3 generally do not exhibit bursts. This 
suggests that the true 15 0(a,7) 19 Ne reaction rate is less than 
that propos ed by CF88. Thi s conclusion is complementary 
to that of Fiske r et all (|2£)06), who were the first to propose 
that the occurrence of type I X-ray bursts is sensitive to the 
strength of the 15 0(a,7) 19 Ne reaction rate, and who found 
that the existence of bursts in systems with / acc « 0.1 suggests 
a lower bound on this rate. This lower bound corresponds to 
/rp w 0.05 in our notation. We also note that a relatively low 
15 0(a,7) 19 Ne reaction rate is consistent with the existence 
of carbon-triggered superbursts, since a low rate would in- 
crease the carbon y i eld resulting from stable nuclear b urning 
ISchatz et alJl!999t ICooper etal1l2006t iFisker et alJ 12006). 
The systems GX 17+2 and Cyg X-2 are exceptions to this 
empirical rule, however, for they e xhibit bursts at accretion 
rates near t he Eddin gton Umit (e.g..|Ka hn & Grindlavl ll984t 
iTawara et all 19841 iKuulkers et alJl2002l) . However, it is pos- 
sible that they show bursts for other reasons, e .g., they har- 
bor m ass donor stars with H-deficient envelopes ( Coop er et alJ 
2006). 

Although the resonant contribution of the 4.03 MeV state 
dominates the 15 0(a,7) 19 Ne reaction rate for the temperatures 
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T < 5 x 10 K at which helium burning can possibly trig- 
ger a type I X-ray burst, other resonances such as the 4.38 
MeV state contribute as well. However, these resonances 
contribute much less than a tenth of the 4.03 MeV contribu- 
tion (e.g., Davids et al 2003), whereas a reduction of the total 
15 0(a,7) 19 Ne reaction rate by a factor of only w 10 is needed 
to stabilize nuclear burning. Therefore, it is the 4.03 MeV a 
width alone that determines the stability of nuclear burning at 
high accretion rates. 

5. CONCLUSIONS 

Using the g lo bal linear stability analysis of 
iCooper & Naravaa d2005h . which is an expanded and 
improved version of the model of NH03, we have investi- 
gated the effects of the hot CNO cycle breakout reaction 
15 0(a,7) 19 Ne on the occurrence of type I X-ray bursts on 
accreting neutron stars at high accretion rates. For low values 
of ftp < 0.1, where f^, is a dimensionless factor by which 
we multiply the ls O(a,7) 19 Ne reaction rate of CF88, the 
hot CNO cycle breakout reaction slightly lowers the critical 
accretion rate above which delayed mixed bursts occur, but 
otherwise the breakout reaction has little effect on the burst 
onset. The predictions of these models are in good agreement 
with observations. For / lp > 0.1, a new regime of prompt 
mixed bursts appears at accretion rates above the rates at 



which delayed mixed bursts occur, and the bursts survive 
up to roughly the Eddington limit. The existence of this 
prompt mixed burst regime up to the Eddington limit is 
consistent with nearly all previous theoretical models. Our 
results support the hypothesis of CN06 that the discrepancies 
between the results of the global linear stability analysis of 
NH03 and the results of time-dependent multi-zone models 
with large reaction networks may be caused by the latter 
models assuming too large a strength for the 15 0(a,7) 19 Ne 
reaction rate at temperatures T < 6 x 10 8 K. The fact that 
observations agree much better with the results of NH03 and 
CN06 implies that the true reaction rate is lower than the 
rate assumed in these multi-zone models. Specifically, we 
suggest that the true a width r Q of the 4.0 3 MeV state in 
19 Ne i s lower than the a width proposed by Langan ke et alJ 
( 1986). Calculati ons us ing multi-z one models such as those 
of lWooslev et alJ fe004l) and lFisker et alJ J2006I) . but with the 
15 0(a,7) 19 Ne reaction rate lowered in strength, need to be 
carried out to either confirm or refute this suggestion. 
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